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By applying the standard adiabatic approximation and using the accurate analytical expression for the 
corresponding local chemical potential obtained in our previous work [Phys. Rev. A 75, 063610 (2007)] 
we derive an effective ID equation that governs the axial dynamics of mean-field cigar-shaped condensates 
with repulsive interatomic interactions, accounting accurately for the contribution from the transverse de- 
grees of freedom. This equation, which is more simple than previous proposals, is also more accurate. 
Moreover, it allows treating condensates containing an axisymmetric vortex with no additional cost. Our 
effective equation also has the correct limit in both the quasi- ID mean-field regime and the Thomas-Fermi 
regime and permits one to derive fully analytical expressions for ground-state properties such as the chem- 
ical potential, axial length, axial density profile, and local sound velocity. These analytical expressions 
remain valid and accurate in between the above two extreme regimes. Following the same procedure 
we also derive an effective 2D equation that governs the transverse dynamics of mean-field disk-shaped 
condensates. This equation, which also has the correct limit in both the quasi-2D and the Thomas-Fermi 
regime, is again more simple and accurate than previous proposals. We have checked the validity of our 
equations by numerically solving the full 3D Gross-Pitaevskii equation. 
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; I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sates (BECs) of dilute atomic gases confined in optical and 
magnetic traps Ol 0, [sll has opened new opportunities for 
investigating the coherence properties of degenerate quan- 
tum systems. From a theoretical point of view, under the 
usual experimental conditions these systems can be accu- 
rately described by the Gross-Pitaevskii equation (GPE) 
10], a mean-field equation of motion governing the behav- 
ior of the condensate wave function '4){t, t) 

\ ih^=[-^V^ + V{v) + gN\i,\'^^l,. (1) 

In the above equation N is the number of atoms, g = 
ATih^a/m is the interaction strength, a is the s-wave scat- 
tering length, and V^(r) is the potential of the confining 
trap. In what follows we shall restrict ourselves to the usual 

\ case of condensates with repulsive interatomic interactions 

■ (a > 0). 

" The GPE has proved to be very successful in describing 
the evolution in time of dilute quantum gases near the zero- 
temperature limit. From a mathematical point of view this 
equation is a time-dependent nonlinear differential equa- 
tion. Since no explicit analytical solutions are known, in 
general, Eq. ([T]l has to be solved numerically. This is a non- 
trivial numerical task that demands a considerable compu- 
tational effort. Moreover, in many circumstances the super- 
fluid dynamics of a zero-temperature condensate can be- 
come chaotic which requires large basis or grid point sets to 
guarantee convergence |5]. In recent years there has been 
particular interest in BECs confined in highly anisotropic 
traps flSBISlllOllilllll. In such geometries the con- 
densate is so tightly confined in the radial or the axial di- 



mension that the corresponding dynamics becomes effec- 
tively one dimensional or two dimensional, respectively. 
The time evolution of these systems with reduced dimen- 
sionality is characterized by two very different time scales. 
Even though usually one is only interested in the evolution 
of the slow degrees of freedom in the effective mean field 
induced by the fast degrees of freedom, from a computa- 
tional point of view it is also required to resolve accurately 
the irrelevant fast degrees of freedom. It is clear that for 
sufficiently anisotropic traps this can represent a computa- 
tional challenge. It is therefore convenient to develop theo- 
retical models that permit one to study the condensate dy- 
namics in terms of effective equations of lower dimension- 
ality. In this regard various approaches have been followed 
in recent years filQIilllEEIill- Among them, the 
effective ID and 2D nonpolynomial nonlinear Schrodinger 
equations by Salasnich et al. ifisll have proved the most 
efficient. 



In this work we derive effective ID and 2D wave equa- 
tions that govern the dynamics of mean-field cigar-shaped 
and disk-shaped condensates with repulsive interatomic in- 
teractions. These equations which incorporate properly the 
contribution from the fast degrees of freedom have the cor- 
rect limits in both the TF and the perturbative regime. Even 
though the equations by Salasnich et al. are in general very 
accurate, we demonstrate that our effective equations are 
more accurate. Moreover, as a consequence of its simplic- 
ity, they also permit one to obtain fully analytical expres- 
sions for various relevant ground-state properties such as 
chemical potentials, condensate lengths, density profiles, 
and local sound velocities. 
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II. CIGAR-SHAPED CONDENSATES 



where we have defined 



Consider a BEC confined in a highly elongated trap. The 
high anisotropy of the trap has important consequences 
on the condensate dynamics which under usual conditions 
becomes governed by two very different time scales. In 
these circumstances the characteristic evolution time of the 
fast transverse motion (~ wj^) is so small in compari- 
son with the characteristic time scale of the axial motion 
(~ w^^) that one can assume that at every instant of time 
the transverse degrees of freedom adjust instantaneously to 
the lowest-energy configuration compatible with the axial 
configuration occurring at that time (adiabatic approxima- 
tion). This implies, in particular, that the correlations be- 
tween transverse and axial motions can be neglected and 
the condensate wave function can befactorized as fBLfll 



= ip{rjL;ni{z,t))(l){z,t), 



(2) 



where = {x,y) and 77-1(2;, t) is the local condensate 
density per unit length characterizing the axial configura- 
tion 



ni{z,t) = N 



\ij{r^,z,t)\' 



Normalizing the transverse wave function to unity 

d^r_L|(/?(ri;ni)|2 = 1, 

Eq. (O takes the desirable form 

n,{z,t) = N\()){z,t)\''. 



(3) 



(4) 



(5) 



Substituting now the wave function ^ into the GPE and 
assuming the confining potential to be separable as V{r) = 
V±{r±^) + Vz{z), one obtains 



V^[z)4) 93(r_L;ni) 



dt 2m dz^ 



4>{z,t). 



(6) 



Both the axial and time variations induced in the transverse 
wave function ip by the axial density rii have been ne- 
glected in the above equation. Neglecting the time deriva- 
tive of if is in fact the essence of the adiabatic approxima- 
tion already discussed. Neglecting the second derivative of 
ip with respect to z requires the axial density to vary suf- 
ficiently slowly along the axial direction. For condensates 
in highly elongated traps such condition usually holds in 
most cases of practical interest. 

Multiplying Eq. ^ by (/?*(r^;ni) and integrating on 
the transverse coordinates r 1 we arrive at 



ih 



dt 



2m dz'^ 



+ V,{z)^ + H±{ni)(j), (7) 



(8) 

Equation (|7]l is an effective ID mean-field equation that 
governs the axial dynamics of cigar-shaped BECs, incorpo- 
rating the contribution from the transverse degrees of free- 
dom through /i^(ni). Substitution of Eq. (|7]l into Eq. ^ 
yields 



^']_ + Vj_{rjL) + gni\ip\ ] ip = fi^{ni)ip 



2m 

(9) 

which shows that at every instant of time and z plane the 
transverse wave function ip satisfies the stationary GPE 
of an axially homogeneous condensate characterized by a 
density per unit length ni{z, t), with ^^{rii) being the cor- 
responding (transverse) local chemical potential. Clearly, 
the usefulness of the effective ID axial equation ^ de- 
pends, to a great extent, on the possibility of finding a sim- 
ple way of solving Eq. (|9l. In this regard, it is important 
to note that rii enters the equation above as a mere external 
parameter, so that, in practice, the solution of the trans- 
verse equation does not require the knowledge of the axial 
evolution. Before considering this problem in more detail 
we will first take a look at the condensate stationary states, 
which must satisfy 



Substituting in Eq. ^ one obtains 



(10) 



2m dz-^ 



+ V;(z)0o + ^±(ni)(/>o = /i(/>o, (11) 



where the condensate chemical potential ^ is the La- 
grange parameter that guarantees the normalization condi- 
tion J dz\(f){z, t)p = 1. When (j)^) varies so slowly that the 
typical length scale of its spatial variations along z is 
much greater than the corresponding axial healing length, 
i.e.. 



A, > 



(12) 



then the first term on the left in Eq. ( fTTl) can be neglected 
in comparison with the mean-field interaction energy and 
one arrives at the local density approximation 



/X = /x^(ni) + V^{z). 



(13) 



When this equation is applicable, the knowledge of the 
transverse chemical potential permits one to obtain an ana- 
lytic expression for the ground-state axial condensate pro- 
file ni{z). In any case, as Eq. (|7]l shows, the analytic deter- 
mination of /x^(ni) is the key ingredient to derive a useful 
ID effective equation of motion for the condensate axial 
dynamics. 
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We shall concentrate in what follows on cigar-shaped 
condensates confined in the radial direction by an axisym- 
metric harmonic potential characterized by an oscillator 
length a± = JUjmuo^, 



y(r) 



(14) 



Introducing dimensionless variables = r^ja^ and 
Tp = a^(p, the transverse equation ^ determining the lo- 
cal chemical potential Ji^ = / tiw±_ takes the form ^ 



]^ 2 1 

--V^ + -r\ + inarii \Tp\^ ] Tp = Jij^{ni)Tp. (15) 



This equation, that depends on the sole parameter arii, 
can be analytically solved in two limiting cases. When 
aril ^ 1 the mean-field interaction energy can be treated 
as a weak perturbation. In this perturbative regime the con- 
densate wave function that minimizes the energy functional 
is given, to the lowest order, by the Gaussian ground state 
of the harmonic oscillator and the condensate is tightly 
confined in the radial direction. Under these conditions the 
radial motion is frozen out, restricted to zero-point oscilla- 
tions, and the condensate becomes effectively one dimen- 
sional. Thus the perturbative regime corresponds to quasi- 
ID mean-field condensates with a local chemical potential 
given by 

7I^(ni) = 1 + 2ani. (16) 

When aril ^ 1 the kinetic energy term can be safely ne- 
glected in comparison with the mean-field interaction en- 
ergy. This is the Thomas-Fermi regime, in which many 
modes of the transverse harmonic trap become excited and 
the condensate exhibits a parabolic radial profile 



1 



Airarii 



1 



(17) 



where the local chemical potential ensuring normalization 
is now given by 



1. (18) 

Note that in passing from the Thomas-Fermi regime to the 
perturbative regime the system undergoes an effective di- 
mensional crossover from a 3D cigar-shaped condensate to 
a quasi- ID mean-field condensate. 

In previous works [20], by using a suitable approxi- 
mation scheme, we derived general approximate formu- 
las that provide with remarkable accuracy (typically bet- 
ter than 1%) the ground-state properties of any mean-field 
scalar Bose-Einstein condensate with short-range repulsive 
interatomic interactions, confined in arbitrary cylindrically 
symmetric harmonic traps, and even containing a multi- 
ply quantized axisymmetric vortex. This approximation 
scheme essentially represents an extension of the Thomas- 
Fermi approximation that incorporates conveniently the 
zero-point energy contribution. In the limiting cases of 



cigar-shaped and disk-shaped condensates the ground-state 
properties follow from explicit analytical formulas that re- 
duce to the correct analytical expressions in both the TF 
and the perturbative regimes, and remain valid and accurate 
in between these two limiting cases, thus accounting prop- 
erly for the corresponding dimensional crossover. Physi- 
cal quantities such as the condensate radius, axial length, 
chemical potential, mean-field interaction energy, kinetic 
and potential energies, density profiles, and local sound ve- 
locities can be easily and accurately obtained in this way 
in terms of the physically relevant parameters (number of 
atoms, trap aspect ratio, and vortex charge). In particular, 
in Ref. [20], we found that the transverse local chemical 
potential of a condensate with an axisymmetric vortex of 
charge q, as a function of the condensate density per unit 
length ni, is given by 



with 



/9, 



22l'Jl(|g|!)2 
(2kl)! 



(19) 



(20) 



The parameter /3~^ accounts for the dilution effect that 
the centrifugal force associated with the vortex has on the 
condensate mean density [|20ll . In the absence of vortices 
g = — > /3q = 1, the equation above simplifies to 

71^ (ni) = Vl + 4ani. (21) 

Clearly, in the appropriate limits, this equation reduces to 
the quasi- ID and TF expressions ( fT6l) and ( fTS]) , respec- 
tively. As we shall see, it also describes correctly the cor- 
responding dimensional crossover. 

The above equations were derived in Ref. ll20ll by using a 
suitable TF-like ansatz for the local density \'\l)(r z)^ of 
a harmonically trapped BEC in its ground state. Such local 
density, which is not factorizable in general, is defined in 
a volume that corresponds to the usual TF ellipsoidal den- 
sity cloud conveniently truncated in order to account for the 
zero-point energy contribution. Note that even though the 
expression ([T9l l was originally obtained for a condensate in 
its ground state (compatible with an axisymmetric vortex 
of charge q) and axially confined by a harmonic potential, 
it is, however, of general validity. This is so because, as al- 
ready said, the solution of the transverse equation ([Tsl l does 
not depend on the particular axial evolution. To convince 
the reader that this is the case we have numerically solved 
Eq. ( ITSl ) for a wave function of the form 

^(rj_) = exp(zg6')^q(r^), (22) 

with q = 0, 1, 2, and 4. The results of the numerical cal- 
culation are shown in Fig. [T] (open circles) along with the 
theoretical prediction obtained from Eq. ( fT9l ) (solid lines). 
As is apparent, Eq. ([T9l l accurately accounts for the dimen- 
sional crossover, the agreement with the numerical results 
being excellent for any value of the dimensionless interac- 
tion parameter aui. The maximum error is smaller than 
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FIG. 1: (Color online) Theoretical prediction for the transverse 
local chemical potential JIj^ as a function of ani for different 
vortex charges g (solid lines). The open circles are exact numer- 
ical results obtained by solving the stationary Gross-Pitaevskii 
equation ( I15l l. 



1% for q = and 2, and smaller than 1.2% for q = 1. 
Even though one expects [20] this maximum error to in- 
crease with q, it is still smaller than 2.5% for q = A. This 
demonstrates that the local chemical potential as given by 
Eq. ( [T9l l is an accurate solution of the transverse equa- 
tion (|9|. This result is of interest in its own right. For 
instance, in Ref. lll9ll . it has been shown that the velocity 
of sound in an axially uniform Bose-Einstein condensate 
immersed in a ID optical lattice and radially confined by 
a harmonic trap can be obtained, in a wide range of opti- 
cal lattice depths, from the knowledge of the local chemi- 
cal potential Jl^ as a function of the linear average density 
rii = N/ d, with d being the lattice period. In such an ap- 
proach the effect of the lattice is incorporated into a suitable 
renormalization of the mass and of the coupling constant. 
In the present work, however, the main interest of Eq. ( [T9l l 
comes from the fact that when substituted in Eq. (Q it leads 
to the following effective ID mean-field equation: 



ih 



dt 



2m dz^ 



+ V,{z)(t> + hw^^[3l + AaN\ 



(23) 

where (3q is given by Eq. ( |20l l and we have absorbed the 
constant term proportional to (|g| + 1) — /3g into the defini- 
tion of the generic axial potential Vz{z). For (7 = 0,1, 2, ... 
the parameter jSq takes the values 1,2,8/3,..., respec- 
tively. 

Equation (|23] | is the equation we were looking for. It 
governs the axial dynamics of arbitrary mean-field cigar- 
shaped condensates with repulsive interatomic interactions 
even in the presence of an axisymmetric vortex of charge 
q, accounting for the effects from the transverse degrees 



of freedom through the term proportional to hio^. Since 
vortices with q > 2 axe dynamically unstable and decay 
into an array of singly quantized vortices iH |2l]], in such 
cases the applicability of Eq. ( |23l l is restricted to times 
shorter than the corresponding decay time. As is apparent, 
the equation above is much simpler than previous propos- 
als. It is also more accurate, as we shall see. Since we have 
found an expression for the transverse local chemical po- 
tential that is very accurate, the validity of Eq. ( |23l > relies 
almost exclusively on the validity of the well established 
adiabatic approximation. 
When Aarii <^ 01 one enters the quasi- ID mean-field 



regime. In this case, Eq. 



with 



2m dz 



5iD = /3, ^ 



reduces to 



+ V,{z)(t> + g^^N\(t>Y ct>, (24) 



27rai 



(3-^2ahuJi 



(25) 



This equation generalizes the well-known ID GPE to the 
case of condensates containing an axisymmetric vortex. In 
the absence of vortices, f3q = \ and one recovers the usual 
equation. When the number of particles is high enough that 
Aarii ^ (3"^, the condensate enters the TF regime. In this 
case, Eq. (E3] | reduces to 



ih 



dt 



2m dz'^ 



+ V,iz)(l) + 2huj^VaN\(l)\(j), (26) 



which again is the correct result as follows from the direct 
substitution of Eq. (fTSl l into Eq. (|7]l. 

The effective ID equation ( |23]) also permits deriving an 
accurate analytical expression for the stationary axial den- 
sity profile ni{z) = N\4>q{z)\'^. To see this we shall 
consider a harmonic axial confinement Vz(z) = T^mujlz^ 



with corresponding oscillator length = ^Jh/moj^. For 
cigar-shaped condensates the trap aspect ratio A = uj^/uj^^ 
must satisfy the inequality A ^ 1. Under these cir- 
cumstances the condition (fT2] | can be easily fulfilled. It 
would be sufficient (though not necessary) that ^ 
a^/ ^JAanl. Taking into account that, in the stationary 
state, A2 is of the order of the condensate axial length it 
is clear that, except for extremely small values of arii, in 
cigar-shaped condensates this condition always holds. As 
a consequence, the stationary equation ( fTTl) reduces to the 
local density approximation ( fT3]) . Substitution of Eq. ( [T9l l 
into Eq. ([Tsl l yields 



AanAz) 



i\q\ + l)+P, 



■q-liVX-zfj 



(27) 



where z = z /a^ . The dimensionless axial half-length Z 
Z/uz follows from the condition ni{Z) = 



^ = (M+l) + i(VAZ,'. 



(28) 
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Substituting this expression in Eq. ( l27l) we obtain 



ni{z) = (3, 



4a V 16a V 

(29) 

with ni(z) = for |z| > Z. In order for this equation to 
be useful one also needs an analytic expression for Z or, 
equivalently, for \i. From the normalization condition 

N= dzni{z), (30) 
J-z 

one finds that the axial half-length satisfies the quintic 
polynomial equation 



Xi, 



(31) 



where Xi = XNa/a± is (apart from the vortex charge) the 
only relevant parameter. An approximate solution of the 
above equation is given by 



Vxz 



1 



+ 



1 



+ 



1 



(15;^j5 + i 57x1 + 345 (3xi//3, 



(32) 

This expression satisfies Eq. (|3T] | for any Xi £ [0, oo), 
with a residual error ll2Cill that is smaller than 0.75% for 
g = and smaller than 3.2% for 1 < |g| < 10. The (lo- 
cal) axial (first) sound velocity Cid of a cigar-shaped con- 
densate, is defined by 

2 ni dnj_ 



'ID 



m drii 

Substituting Eq. ( [T9l l in Eq. ( |33l l one obtains 



mc 



ID 



I 4a'^nl{z) 
(3l + ^an,[z) 



(33) 



(34) 



Equations (I28l) - (l34l) c oincide with the results we obtained 
in previous works ll20ll . Note, however, that the derivation 
is quite different. The formulation of Ref. llOh is not re- 
stricted to condensates of specific geometry. The formulas 
obtained there are valid for condensates confined in arbi- 
trary cylindrically symmetric harmonic traps and reduce 
to those derived in the present work in the limiting case 
of highly elongated condensates. In particular, Eq. ( |28l l 
above appears in Ref. llioll as part of the initial ansatz for 
the local density \^{r^, z)\^ . On the other hand, Eq. ( |29l l 
follows in Ref. L2Q1 from a direct integration of this ansatz 
over the radial coordinate, while Eq. dSTJ is obtained after 
integrating \'ip{r^, z)\'^ over both the radial and the axial 
coordinates. 

As shown in Ref. llioll . Eqs (l28ll-(l34ll predict very ac- 
curately the condensate ground-state properties. They also 
reduce to the correct expressions in the two analytically 
solvable regimes. In particular, in the quasi- ID mean-field 
regime (corresponding to Xi *C 1 ^ clui <C 1) the axial 
density profile ( |29l l is well approximated by the first term 



on the right-hand side while, in the TF regime (correspond- 
ing to Xi 1 — ^ CL^i ^ 1), the last term on the right- 
hand side is the only one that contributes significantly, in 
good agreement with previous results obtained in these two 
particular limits [12]. Equation ( |29l l also reproduces accu- 
rately the axial density profile in between these two limit- 
ing cases. This has been demonstrated in Ref. [20], where 
we compare the theoretical prediction obtained from Eqs. 
(l29l l and (l32l l with exact numerical results obtained from 
the full 3D GPE. The agreement is always very good except 
at the condensate edges, as expected for a TF-like expres- 
sion obtained by neglecting the derivative kinetic term. In 
this regard, we note again that the derivation of the above 
equations followed in the present work is based on the lo- 
cal density approximation which in turn requires the con- 
dition (fT2] l to be true. As already said, however, for cigar- 
shaped condensates with A <C 1 this requirement is not 
very demanding and can be satisfied even in the perturba- 
tive regime (see also Ref. (v^). 

Our effective ID equation (l23l l is a Schrodinger equation 
with an effective mean-field potential 



(35) 



As mentioned in the Introduction, by using a variational ap- 
proach ['22!], Salasnich et al. obtained in Ref. Iil5i1 effective 
ID and 2D nonpolynomial nonlinear Schrodinger equa- 
tions (NPSEs) for the axial and radial dynamics of cigar- 
shaped and disk-shaped condensates, respectively. Their 
ID equation coincides with our equation above but with 
a different effective mean-field potential 




Next we will compare our effective ID equation ( 1231) with 
that by Salasnich et al. These authors demonstrated that 
their equations provide much more accurate results than 
any other proposed effective equation. Actually, their re- 
sults are always very close to the exact results obtained 
from the 3D GPE. It is thus sufficient to compare our 
equations with those proposed in Ref. ifTsll . As we shall 
see, our effective equation is more accurate. Moreover, 
it allows treating condensates containing an axisymmet- 
ric vortex, with no additional cost. For instance, to ac- 
count for a q = 1 vortex it suffices to make the change 
(3q=Q = 1 — > /3q=i = 2 in Eq. (l35l l. Our effective equation 
also permits one to derive useful and accurate analytical 
expressions for ground-state properties such as the chem- 
ical potential, axial length, axial density profile, and local 
sound velocity. In this regard, it should be noticed that 
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FIG. 2: (Color online) Axial density profile ni (z) of cigar-shaped 
condensates in a harmonic trap with aspect ratio A — ujz/ijjx^ = 
0.1 and xi ~ 0.1, 1 and 10. (a)-(c): Ground-state {q — 0) equi- 
librium configuration, (d): Equilibrium configuration compatible 
with a vortex of charge q — 1 and 2. Solid lines are numerical re- 
sults obtained from our effective ID equation l |23l l. Dashed lines 
are numerical results obtained from the ID NPSE by Salasnich et 
al. Open symbols are exact numerical results obtained from the 
full 3D GPE. 
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FIG. 3: (Color online) Time evolution of the mean squared am- 
plitude (z^) of cigar-shaped condensates with Xi — 1 ^nd 10 and 
vortex charges g = and 1 after a perturbation cjz — > I-Il^Jz- 
The aspect ratio of the harmonic trap before perturbation is 
A = LJz/i-^± = 0.1. Solid lines are numerical results obtained 
from our effective ID equation | |23I >. Dashed lines are numeri- 
cal results obtained from the ID NPSE by Salasnich et al. Open 
symbols are exact numerical results obtained from the full 3D 
GPE. 



although Salasnich et al. also obtained in Ref. ll23h an- 
alytical expressions for the axial density profile and local 
sound velocity, such expressions depend on the condensate 
chemical potential which because of the complexity of Eq. 
(l36l l in general cannot be determined analytically except in 
the two analytically solvable regimes or in the simple case 
of a homogeneous condensate with no axial confinement. 
Finally, note that the approach of Ref. [[15] is also valid for 
condensates with attractive interatomic interactions. 

Figure |2] shows the axial density profile ni(z) of cigar- 
shaped condensates in a harmonic trap with aspect ratio 
A = uJz/uJi_ = 0.1 in the different relevant regimes, 
which, for a given vorticity q, are completely characterized 
by the value of the sole parameter Xi- Figures I2a)42c) 
correspond to the ground-state (q = 0) equilibrium config- 
uration, while Fig. ^d) shows the equilibrium configura- 
tion compatible with a vortex of charge q = 1 and 2. Solid 
lines are numerical results obtained from our effective ID 
equation (|23] |. Dashed lines are numerical results obtained 
from the ID NPSE by Salasnich et al. and open symbols 
are exact numerical results obtained from the full 3D GPE. 
As is apparent from Fig. |2a), in the perturbative regime 



(Xi = 0.1), both approaches are practically indistinguish- 
able, a consequence of the fact that both have the correct 
perturbative limit. However, as Xi increases our effective 
ID equation is always more accurate [Figs. |2b) and|2c)], 
which is a consequence of the fact that our approach, un- 
like that by Salasnich et al., also reproduces correctly the 
TF limit. The difference, though small, can be clearly ap- 
preciated even for Xi = 1- Moreover, our approach also 
allows treating the case of condensates containing an ax- 
isymmetric vortex, with no additional cost. Figure l^d) 
shows the equilibrium configuration of condensates with 
a vortex of charge q = 1 and 2, obtained from the same 
effective ID equation ( |23l l by simply taking (3q = 2 and 
8/3, respectively. As is apparent, the agreement with the 
exact numerical results (open symbols) is again very good. 

Our approach also reproduces very accurately the time 
evolution of cigar-shaped condensates. To see this we start 
from the equilibrium configuration of condensates confined 
in a harmonic trap with aspect ratio A = u)z/u)±^ = 0.1. 
Then we introduce a sudden perturbation by increasing the 
axial confinement frequency as uj^ ^-^^z and follow 
the subsequent evolution in time of the mean squared axial 
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FIG. 4: (Color online) Evolution of the axial density profile ni (z) 
of a cigar-shaped condensate in a harmonic trap with aspect ratio 
A = LUz/uj± — 0.1 and xi ~ 2 for different times r = cozt. Solid 
lines are numerical results obtained from our effective ID equa- 
tion i23l . Dashed lines are numerical results obtained from the 
ID NPSE by Salasnich et al. Open symbols are exact numerical 
results obtained from the full 3D GPE. 



amplitude (z^) = J d'z'z'^ni{z,t). Figure[3]shows the cor- 
responding results obtained for condensates with Xi = 1 
and 10 and vortex charges g = and 1. As before sohd 
hnes are numerical results obtained from our effective ID 
equation (l23l l. dashed lines are numerical results obtained 
from the ID NPSE by Salasnich et ah, and open symbols 
are exact numerical results obtained from the full 3D GPE. 
For Xi = 0-1 (not shown in the figure) both approaches 
produce results that are indistinguishable on the scale of 
the figure from the exact results. It is apparent that our ap- 
proach is more accurate than that by Salasnich et al. More- 
over, as Fig. [3la) shows, it also produces very accurate 
results for condensates containing an axisymmetric vortex. 
Note also that for Xi = 1 and g = 0, Fig. Oa) shows that 
the results from the ID NPSE begin to dephase at large 
times. To see this behavior more clearly, it is convenient to 
consider a more complex situation. To this end we consider 
a BEC with Xi = 2 in a harmonic trap with A = 0.1, sub- 
ject to a blue-detuned laser beam modelled by the Gaussian 
potential 



(37) 



with Vo = 12^2 and Zq = v2Aa2. We let the con- 
densate reach its equilibrium configuration and then, at 



r = uizt = 0, we suddenly switch off the laser beam and 
let the system evolve in the harmonic trap. Figure |4] shows 
the time evolution of the axial density profile ni(z). Solid 
lines are numerical results obtained from our effective ID 
equation (1231) . dashed lines are numerical results obtained 
from the ID NPSE by Salasnich et ah, and open symbols 
are exact numerical results obtained from the full 3D GPE. 
At r = 1 one recognizes in each half-axis a dark soliton- 
like structure that develops from the initial configuration 
and propagates toward the corresponding condensate edge 
(only the positive half-axis is represented in the figure). At 
r = 2 the soliton becomes black and it comes back to- 
ward the condensate center for r > 2, as can be seen in the 
snapshot at r = 3. This figure demonstrates again that our 
approach is more accurate than that by Salasnich et al. In 
particular, it is apparent that the results obtained from the 
ID NPSE become somewhat dephased with respect to the 
exact results. 



III. DISK-SHAPED CONDENSATES 

Consider now a condensate confined in an anisotropic 
trap that is much stronger in the axial than in the radial 
direction. In this case, using similar arguments as before, 
one can resort to the adiabatic approximation and assume 
that at every instant of time the (fast) axial degrees of free- 
dom adjust instantaneously to the equilibrium configura- 
tion compatible with the radial configuration occurring at 
that time. Under these circumstances the condensate wave 
function can be factorized as 

V'Cr, t) = ip{r^,t)^{z; n2(r^, t)), (38) 

where 71-2 (r^,t) is the local condensate density per unit 
area characterizing the radial configuration 

n2{v^,t)=Njdz\^{v^,z,t)\\ (39) 

The normalization condition 

^ dz\<j){z-n2)\^ = 1 (40) 

leads to 

n2{v^,t) = N\^{v^,t)\\ (41) 
After substituting Eq. ( [38]) into Eq. ([Tl, one arrives at 

dt 2m 

-^-g^ + yz{z)(j) + gn2{r^,t) \<l)\^ (/)] ip{r^,t), 

(42) 

where as before we have assumed the confining potential 
to be separable. Integrating out the fast degrees of freedom 
one obtains 

ih-^ = -—Vlip + V^{r^)ip + Hz{n2)ip, (43) 



if>'-^ + 77— Vj.'^ - V^{r^)(p ] (t>{z;n2) 
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with 



dz ( 



+ V,{z)+gn2\(l)f 



2m dz'^ 



Substituting Eq. (|43] | in Eq. (|42]| one finds 



(46) 



(44) 



2m dz^ + ^"^^^ ^ ^""^ J " ^-("2)<A- (45) 

The effective 2D mean-field equation (|43l l governs the ra- 
dial dynamics of disk-shaped BECs, incorporating the con- 
tribution from the axial degrees of freedom through the (ax- 
ial) local chemical potential fj,z{n2). This quantity, in turn, 
follows from the stationary Gross-Pitaevskii equation i 
determining the axial configuration (p{z; 712). 
The condensate stationary states 

satisfy 

- 7;— V^(/7o + V^{r^)ipo + fi,{n2)ipa = fi(po, (47) 
2m 

where the chemical potential fi guarantees the condition 
/ (Fr±\(p{r±,t)\'^ = 1. When the typical length scale 
of the spatial variations of (po is much greater than the cor- 
responding radial healing length, i.e., 

> ^ (48) 

yj2m{n^{n2) - ^hwz) 

then Eq. (l47l) reduces to the local density approximation 

^ = ^.(n2) + 14(rx). (49) 

We shall consider in what follows disk-shaped condensates 
confined in the axial direction by a harmonic potential char- 
acterized by an oscillator length = y^h/mtUz, 

V{r) = ^mu;y + V^{r^). (50) 

In terms of dimensionless variables 1 = z /oz, cp = ^/oT'^' 
and Jl^ = ^z/hwz, the axial equation (l45l) reads 



1 5' 1-2 . ITI- 

"2^ + 2^ +4vraa,n2|(/)| 



cf) = //z(n2)(/'. 



(51) 

In Ref. ilOll . we found the following expression for the 
axial local chemical potential 'p.^ as a function of the con- 
densate density per unit area 712 : 

1 n 2 



lJ'An2) 



where ry 
with 



4 + 6^1 — + 247raazn2 and ^1 



(52) 
{^2 - 1) 



1^2 J 



E V2A + e(n2 -0.1) 

1 - ^) (1 - (10n2)-'/') . 




z 2 



FIG. 5: (Color online) Theoretical prediction for the axial local 
chemical potential Jt^ as a function of aazn2 (solid lines). The 
open circles are exact numerical results obtained by solving the 
stationary Gross-Pitaevskii equation dSll l. 



In the above equation, G(x) is the step function and n2 = 
002722 is the only relevant parameter. Note that in Ref. 
[20] we used the same expression i53[ but in terms of the 
dimensionless parameter X2 = A^a/A^Oz instead of n2. 
It was shown there that X2 is the only relevant parame- 
ter for the description of harmonically trapped disk-shaped 
condensates, so that n2 and X2 are directly related to each 
other. In particular, the quasi-2D perturbative regime cor- 
responds to X2 1 ^ ^2 1; while in the TF regime 
X2 ^ 1 ^2 ^ 1- On the other hand, as explained 
in Ref. flB], the (slowly varying) last term in Eq. (l53l l 
was conveniently introduced to ensure that in the TF limit 
— > 1, its specific functional form being not very rel- 
evant. This way one guarantees the correct expression for 
71^(712) in both the TF and the perturbative regimes. On 
these grounds and taking into account that X2 ~ ^2: we 
have simply made the substitution X2 ^2 in Eq. (l53l l 
which now becomes a very slowly varying function of n2- 

Figure|5]shows the theoretical prediction for /i^ obtained 
from Eq. (|52] | (solid lines) along with the exact results ob- 
tained from the numerical solution of the axial equation 
dST]) (open circles). As is apparent, not only does Eq. (l52l l 
have the correct limit in the two relevant extreme regimes, 
but it also accounts accurately for the corresponding di- 
mensional crossover. 



(53) 



Substituting Eqs. and with n2 = N\(p{rj_,t)\'^ 
into Eq. (l43l l, we finally obtain the desired effective 2D 
equation. This equation governs the transverse dynamics of 
arbitrary mean-field disk-shaped condensates with repul- 
sive interatomic interactions, accounting properly for the 
effects from the axial degrees of freedom. In the quasi-2D 



FIG. 6: (Color online) Radial density profile n2{rj_) of disk- 
shaped condensates in a harmonic trap with aspect ratio A = 
LUz/uj± — 10 and X2 ~ 0.1, 1, and 10. (a)-(c): Ground-state 
(q = 0) equilibrium configuration, (d): Equilibrium configura- 
tion compatible with a vortex of charge q = 1. Solid lines are 
numerical results obtained from our effective 2D equation ( 143 b . 
Dashed lines are numerical results obtained from the 2D NPSE 
by Salasnich et al. Open symbols are exact numerical results ob- 
tained from the full 3D GPE. 



FIG. 7: (Color online) Time evolution of the mean squared ampli- 
tude {r\) of disk-shaped condensates with X2 = 1 and 10 after a 
perturbation uj^ I.Iljj^. The aspect ratio of the harmonic trap 
before perturbation is A = Uyjuj^^ = 10. Solid lines are numeri- 
cal results obtained from our effective 2D equation ( 143 b . Dashed 
lines are numerical results obtained from the 2D NPSE by Salas- 
nich et al. Open symbols are exact numerical results obtained 
from the full 3D GPE. 



mean-field regime (aazn2 <C 1), it reduces to 



ih 



dip 
dt 



2m 



Vlp + V^ir^)ip + g2uN\^\^^, (54) 



where = g/v27raz = 2\/2TraazhuJz and we have 
absorbed a constant huJz/2 into the definition of the generic 
transverse potential V±{rj_). Similarly, in the TF regime 
(aa^n2 ^ 1), the effective 2D equation reduces to 



ih 



dip 
'dt 



2m 



+ hiO, 



{?,TT/^j2)aa,N\ipY 



2/3 



ip. 



(55) 



It can be easily verified that both Eq. ( l54l) and Eq. 
are the correct limits of the underlying 3D GPE. As in the 
ID case, our effective 2D equation permits one to derive 
analytical expressions for ground-state properties of disk- 
shaped condensates such as the chemical potential, con- 
densate radius, radial density profile, and local sound ve- 
locity. These analytical expressions reduce to the correct 
formulas in both the TF and the perturbative regimes, and 
remain valid and accurate in between these two limiting 
cases. Since the calculations are more complicated than 



before we refer the reader to our previous work for details 



Next we shall consider disk-shaped condensates in a 
confining potential that is also harmonic in the radial direc- 
tion = ^mLv\r\. The corresponding oscillator 

length is a± = \/h/muJ±^ and the trap aspect ratio now 
satisfies the inequality A = lo^I'^l ^ 1- In Fig. |6]we 
show the radial density profile n2(r^) of condensates in a 
trap with A = 10, in the different relevant regimes char- 
acterized by the parameter X2- The dimensionless variable 
is defined as = r^j a^. Figures |6ta)-[6tc) corre- 
spond to the ground-state (q = 0) equilibrium configura- 
tion, while Fig. [SJ^d) corresponds to the equilibrium con- 
figuration compatible with a g = 1 vortex. Solid lines are 
numerical results obtained from our effective 2D equation 
( l43l) with /izdv^n given by Eqs. (|52l) and (gSll. Dashed 
lines are numerical results obtained from the 2D NPSE by 
Salasnich et al. ifisll and open symbols are exact numerical 
results obtained from the full 3D GPE. In the perturbative 
regime, for Xi = O-l- the 2D NPSE is somewhat more ac- 
curate. This is precisely the parameter region where the er- 
ror of our effective equation is maximum. As X2 increases, 
our equation becomes more accurate. This can be appreci- 
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ated in Figs. ISb)-[6td) which show that our results remain 
very accurate as X2 increases. On the contrary, the error of 
the 2D NPSE increases with X2> being of the order of 5% 
for X2 = 1 and of the order of 10% for X2 = 10. This 
is a consequence of the fact that the underlying (Gaussian) 
variational wave function that leads to the 2D NPSE cannot 
reproduce properly the TF regime. 

To investigate the validity of our effective equation in 
time -dependent problems we have followed the same pro- 
cedure as before. We start from the equilibrium config- 
uration of condensates confined in a harmonic trap with 
A = 10. We then perturb the system by changing instanta- 
neously the radial confinement frequency as uj± — > 1.1uj± 
and follow the subsequent evolution in time of the mean 
squared radial amplitude (r^) = / (Pr±r\n2{r±,t). 
Figure |7] shows the corresponding results for condensates 
with X2 = 1 and 10. Solid lines are numerical results ob- 
tained from our effective 2D equation (l43l l with Hz{\ip\'^) 
given by Eqs. ( [52l ) and ( [53] l, dashed lines are numerical 
results obtained from the 2D NPSE, and open symbols are 
exact numerical results obtained from the full 3D GPE. As 
is apparent from the figure, our effective 2D equation is 
again more accurate than the 2D NPSE by Salasnich et al. 

IV. CONCLUSION 

In this work, by applying the standard adiabatic ap- 
proximation and using the analytical expression for the 
transverse local chemical potential obtained in our previ- 
ous work 120|], we have derived an effective ID equation 
that governs the axial dynamics of mean-field cigar-shaped 
condensates with repulsive interatomic interactions. This 
equation, which incorporates accurately the contribution 
from the transverse degrees of freedom, has the correct 
limits in both the quasi- ID mean-field regime and the TF 
regime. Since our expression for the local chemical po- 
tential is very accurate, the validity of the above equation 
relies almost exclusively on the validity of the well estab- 
lished adiabatic approximation. 

We have compared our effective ID equation with the 
ID NPSE by Salasnich et al. which provides more accurate 
results than any other previously proposed effective equa- 
tion. We have demonstrated that our effective ID equation 
is more accurate, which, in part, is a consequence of the 
fact that our approach, unlike that by Salasnich et al., re- 
produces correctly the TF limit. Moreover, our equation al- 
lows treating condensates containing an axisymmetric vor- 
tex with no additional cost. Because of its simplicity, it 
also permits one to derive fully analytical expressions for 
ground-state properties such as the chemical potential, ax- 
ial length, axial density profile, and local sound velocity. 
These analytical expressions reduce to the correct analyti- 
cal formulas in both the TF and the perturbative regimes, 
and remain valid and accurate in between these two limit- 
ing cases. 



Following the same procedure we have also derived an 
effective 2D equation that governs the transverse dynamics 
of mean-field disk-shaped condensates, accounting prop- 
erly for the contribution from the axial degrees of freedom. 
This equation is also more accurate than the 2D NPSE by 
Salasnich et al. As in the ID case, from this effective equa- 
tion, which also has the correct limits in both the quasi- 
2D and the TF regime, one can derive analytical expres- 
sions for ground-state properties of cigar-shaped conden- 
sates such as the chemical potential, condensate radius, ra- 
dial density profile, and local sound velocity. 
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